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Abstract 

In this paper we consider quadratic stochastic operators designed 
on finite Abelian groups. It is proved that such operators have the 
property of regularity. 
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1 Quadratic stochastic operators 



Let 



n 



S' 



in— 1 



{x= (x x ,... ,x n )€R n : 



Xi = 1 , Xi > \/i = 1 



• • • ,n} 



(1) 



be the (n — 1) -dimensional simplex in IR n . The transformation V : 




is called a quadratic stochastic operator (q.s.o.), if 
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where 

Pij,k > o, 

Pij,k Pji,k: 

n 

^Pij,k = 1 (3) 
k=l 

for arbitrary i,j,k £ {1, ■ • ■ , re}. Such operators have applications 
in mathematical biology, namely theory of heredity, where the coeffi- 
cients pij t k are interpreted as coefficients of heredity [1-3]. 

Assume {V k x : k = 0, 1, ■ ■ ■ } is the trajectory of the initial point 
x € S™ -1 , where V k+l x = V(V k x) for any k = 0, 1, • • ■ 

Definition 1 A q.s.o. V : S 71 ^ 1 — > S 11 ^ 1 is called ergodic (respectively 

regular) if for any initial point x £ S 1 ™ -1 the limit 

^ k-l 

lim - V l x 

k^oo k ^ 

(respectively the limit lim V k x) exists. 

k^oo 

Evidently, any regular q.s.o. V has the ergodic property, but the 
converse is not necessarily true. 

To determine whether some q.s.o. is ergodic or regular is rather 
complicated problem. 

S. Ulam in [4] presupposed the assumption that any q.s.o. V is 
ergodic. Later M. Zakharevitch [5] showed, that this is an incorrect 
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hypothesis in general. More precisely, he proved that the q.s.o. V, 
which is defined on the simplex 



S 2 = {(x, y,z) : x,y,z>0, x + y + z = 1} 



by the formula 



V : < 



x = x 2 + 2xy; 
y = y 2 + 2yz; 



z = z 2 + 2xz. 

is not an ergodic q.s.o. and respectively is not regular. 

Later in [6] necessary and sufficient conditions were established for 
the ergodicity of the so-called Volterrian q.s.o.: 

x = x{\ + ay — bz); 

^ ' 1 V = 2/(1 — ax + cz); 

z = z(l + bx — cy) 

where a,b,c £ [—1,1]. For a,b,c = 1 we have the initial example of 
Zakharevitch. 

Theorem 1 Any q.s.o. in the above form is non-ergodic if and only 
if the three parameters a, b, c have the same sign. 



Below we'll construct one class of q.s.o. and prove that all such 
q.s.o. are regular and respectively ergodic. 



2 The design of quadratic stochastic 
operators 

Let G be a finite Abelian group and S(G) be a set of all probabilistic 
measures on G. It is evident, that if \G\ = n, then S(G) coincides with 

Let further H C G be a subgroup of G and {g + H : g € G} be 
the cosets of H in G. Assume fi € S'(G) is a fixed positive measure, 
that is fi(g) > for any g G G. Then we define the coefficients Pf g ,h, 
where /, g, h G G in the following way: 



P/9,h = S 



otherwise. 

It is easy to check that for arbitrary f,g,h G G the conditions (3) are 
satisfied. It is also evident that if H = {e}, where e is the neutral 
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element of group G, then 



Pfg,h = S 



1 if h = f + g; 
0, otherwise 



and if H = G, then 



Pfg,h = v(h) Vf,geG 
. In the common case q.s.o. V on S(G) is defined as 

(Vx) h = Pfg,h x f x g 
f,geG 

for all h G G, where x = {xt, t G G} G S(G) and Pf g ,h as above. 
If H = {e} then the q.s.o. V is defined as 

(Vx) h = X f X 9 
f,geG,f+g=h 

and if H = G, the q.s.o. is defined as 

(Vx) h = n(h) 

for arbitrary h G G. 

Let us fix a positive measure fi G <S(G) and subgroup .ff of group 
G. Assume hh is the factor-measure on factorgroup G/H, that is 



hen 



for any g G G and Vh is a q.s.o. on S(G/H), which is defined by 
measure fin- It is easy to show, that the trajectorial behaviour of V 
and Vh are similar, and so it is enough to study the q.s.o. generated 
by the trivial subgroup. 

Below we consider q.s.o. constructed by the trivial subgroup H = 



Let v G S(G) be a Haar measure on G. As G is a finite Abelian 
group, the Haar measure v on G is a uniform distribution on G,that is 
it is a centre of the simplex S(G). We prove the following 

Theorem 2 Almost all orbits tend to the center of the simplex. 

Proof: The following lemma is a key ingredient: 

Lemma 1 HVarHoo < ||x||oo- 

Proof of lemma: Let's define the following function 



Notice that there is no dependence on n. Moreover, / has an explicit 
form: 



{e}. 



3 The main result 



n 




f{p) = k P 2 + (l - kpf if pe[ 



l 



l 



fc + l' k 
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We can prove this formula by induction . 

The sum ^ xj cannot have a maximum in the interior of the do- 
main (by the Lagrange method). So, we can look for the maximum 
at the set x n = p. It follows that 



rt-1 

2 



fn{p)=p 2 + max ^x. 

Let Xi = (1 - p)yi.Then Ylo" 1 Hi = 1 and Hi G [°>i^]- So > tne 
following recurrent formula is true: 



fn(p) =P 2 + (1-P) 2 j n-l(: V 



l-p' 

This allows an easy check of the expected formula. Now it is easy to 
see that 

f(p)-p = kp 2 -p+(l-kp) 2 = (l-kp)(l-kp-p) = (l-kp)(l-(k+l)p) < 0, 
that is 

f(p) < P- 

From Cauchy-Bunyakowski inequality: 



V jeG V j£G j&G 



x) 



If maxxj < p then, f(p) < p implies {Vx)i < p or max(V2;) ! < p. 
So, the norm ||x||oo decreases on the orbits. Lemma is thus proven. 



Actually, the set {x : ||x||oo < p} is mapped to the set {x : ||x||oo < 
f(j))}- We can easily check, that the iterations of the function / tend 
to its fixed points. 

Now consider a point x. We consider the case, when a point in orbit 
does not tend to a centre. It's lo— limit set lies at a set {x : g(x) = j:}. 
This u;— limit set is finite. So, it is nothing else, then a periodic orbit. 



We consider the periodic orbits on the set \\x 




Since HVa^loo = ^, there exists k coordinates of x equal to p 
All other coordinates are zeros. This is also true for Vx. So, there 
exists sets A = {i : Xi = ^} and B = {j : (Vx)j = |}, such that 
|A| = |B| = k 

and 

1. Vi e B \(i - A)f]A\ = k ^ i - A = A, 

2. Vi£ £ A)f]A\ = -f (i-A)n^ = 0- 
From the second item: 

Mi ^ B Vx, y £ A : i — x ^ y; 
Vi (£ B Vj £ A + A, i± j; 
A + AcB. 
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From the first item 

Vi G B 3j, k G A : i = j + k; 
BcA + A. 

So, there exists a set A, such that |A| = |A + A\ = k, 

Vi G A + A: i- A = A. 

Claim: This is equivalent to the following property: 

Vi,j,keA: i + j-k£A. 

Proof: — > Vx, y £ A x + y — A = A. Vx, y, z£ii + |/ - z £ A. 
<— Mx,y G Ax + y- A C A. Vz G i + iz-i C A. Since 
1 2 - A | = z - A = A. 

Va,i,j G A a + i - j G A. G A A + i - j c A. 

e A A + i c A + j. 

So, j e A A + i = A + j, and hence , 

A + A=\J(A + j) = A + i. 

jeA 

So, \A + j4| = k. The claim is proved. 

This property is equivalent to the following : there exists a point 
p G A and a subgroup il, such that A = p + H. So, the points 
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Vx, V 2 x, V 3 x, ■ ■ ■ correspond to the sets 2p + H, 4p + H, 8p + H , • • • . 
Actually, such sequences are pre-periodic. 

We now prove the instability of such periodic orbits. Let I be it's 
period and V 1 = T be the corresponding first return map. Then T 
has an instability direction (an eigenvector, whose eigenvalue is greater 
than 1). It is easy to check, that vector e with coordinates 

1, if s Ei + H; 
- \-Hf sej + H; 
0, other cases 

belongs to the plane {x : ^ Xj = 0} (a tangent plane for the simplex) 
and realizes this instability direction. 

There is no such directions for the centre by the coincidense of 
each two classes. 

So, the basin of attraction for such orbit consists of subvarieties of 
strictly positive codimension.We conclude that almost all orbits tend 
to the centre of the simplex. 

Remark 

We can also consider an infinite-dimensional case. All preceding 
arguments remain true. But there is no natural measure on a sim- 
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plex,and so we cannot use the term "almost all". 

Another disappointment is that "most" of the orbits diverge in 
the l\ topology. It follows from the fact, that for "most" orbits the 
sequence \\V n 

%\ |oo tends to zero. All possible l\ limit points are 



limit points. Since the only limit point is zero, then the subsequence 
of the orbit l\ tends to a point, which does not belong to a simplex. 
This contradiction implies the absense of h limit points. 
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